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We report Hubble Space Telescope observations of 6 gravitational lenses 
^ ■ with the Advanced Camera for Surveys. We measured the flux ratios between 

the lensed images in 6 filters from 8140A to 2200A. In 3 of the systems, 
HE0512-3329, B1600+434, and H1413+117, we were able to construct UV 
extinction curves partially overlapping the 2175A feature and characterize the 
^ ■ properties of the dust relative to the Galaxy and the Magellanic Clouds. In 

cn ■ 

a^ 

lO ' properties of the quasar accretion disk. For a Gaussian model of the disk 

f^ ■ exp(— r^/2r^), scaling with wavelength as r^ oc A^, we estimate rs(A3363) = 4^2 



O 



% 



(7 ± 4) light-days and p = 1.1 ± 0.6 (1.0 ± 0.6) for a logarithmic (linear) prior 
on Tg. The remaining two systems, FBQ0951-I-2635 and SBS1520-I-530, yielded 
no useful estimates of extinction or chromatic microlensing. 



Subject headings: cosmology: observations — gravitational lensing — dust, 
extinction, accretion, accretion disk — galaxies: ISM 



""^ Based on observations made with the NASA/ESA Hubble Space Telescope. The Space Telescope Science 
Institute is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. 
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1. Introduction 



The wavelength-dependent flux ratios between gravitationally lensed images of quasars 
can be used to probe both the extinction law in the lens galaxy and the unresolved structure 
of the quasar. Each image is differentially extincted by dust along the path of the image 
through the lens galaxy, and this can be used to measure the extinction, the extinction law, 
and the lens redshift (e.g. Nadeau et al. 1991; Falco et al. 1999; Motta et al. 2002; Munoz 
et al. 2004; Mediavilla et al. 2005; Ehasdottir et al. 2006, Mosquera et al. 2011). The 
second effect, microlensing by the stars in the lens galaxy (Chang & Refsdal 1979), leads 
to wavelength dependent changes in the flux ratios because the effective size of the quasar 
accretion disk varies with wavelength (Wisotzki et al. 1993, 1995; Wucknitz et al. 2003; 
Anguita et al. 2008; Bate et al. 2008; Eigenbrod et al. 2008; Poindexter et al.2008; Floyd 
et al. 2009; Mosquera et al. 2009, 2011; Blackburne et al. 2010; Mediavilla et al. 2011). 
Detections of "chromaticity" between the images of a lensed quasar are useful for studying 
both phenomena if they can be disentangled. 

In extragalactic astronomy, understanding dust is crucial to understanding galaxies, 
through its effects on estimates of star formation rates and galaxy evolution (e.g. Conroy 
et al. 2009), cosmology, through its effects on SNe la fluxes (e.g. Jha et al. 2006), and the 
interpretation of gamma-ray burst afterglows (e.g. Jakobsson et al. 2004)). Unfortunately, 
classical methods for obtaining accurate extinction curves to characterize dust cannot be 
used outside the Local Group because they depend on detailed measurements of individual 
stars. Gravitational lenses are one of the best quantitative astrophysical probes of dust 
properties at intermediate redshifts given lenses with the right amount of dust and the 
appropriate combinations of redshifts. If there is too little dust, it is difficult to measure the 
extinction at long wavelengths and microlensing is more likely to dominate the chromaticity. 
If there is too much dust, it becomes impossible to measure extinction curves into the 
rest-frame ultraviolet. Similarly, the lens redshift must be high enough to make the 
rest-frame 2175A dust feature observable, while the source redshift must be low enough 
to avoid having the quasar continuum blocked by absorption in the intergalactic medium. 
Similar considerations hold for studying chromatic microlensing over the broadest possible 
wavelength baseline. 

We selected six lenses from the survey of extinction by Falco et al. (1999) that roughly 
satisfied these criteria: HE 0512-3329, FBQ 0951+2635, HE 1104-1805, H 1413+117, 
SBS 1520+530, and B 1600+434. As we report in §2, we observed them in 6 filters spanning 
2200A to 8IOOA (F220W, F250W, F330W, F435W, F555W, F625W and F814W) using the 
Hubble Space Telescope (HST) and the ACS/HRC camera. This approach ensures that we 
can measure the image flux ratios without contamination from the lens or the host galaxy 
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of the quasar. Section 2 also outlines how we model the results to study extinction and 
microlensing. In §3 we present the results, reporting on the extinction curves of three of the 
systems and the chromatic microlensing in one system. Section 4 summarizes the results 
and lessons for future observations. 



2. Observations and Analysis 

Table 1 provides a log of our ACS/HRC observations based on 13 HST orbits in 
Cycle 12, Table 2 summarizes previous HST observations of these systems and Table 3 
presents our new photometry. The observations in each filter consisted of multiple, dithered 
sub-exposures which were corrected for cosmic rays and combined using standard methods. 
We modeled the images following the procedures of Lehar et al. (2000). The images were fit 
as a combination of point sources, de Vaucouleurs and exponential disk profiles convolved 
with TinyTim (Krist & Hook 1997) PSF models. We determined the relative astrometry of 
the components and the structure of the lens galaxy using the CASTLES H-band images 
where the lens galaxy is best detected and characterized (see Fig. 1). These were then held 
fixed and the remaining images were fit to determine the fluxes of the components in each 
filter. For the bluer filters the lens galaxy was undetected and we could easily confirm the 
model fits with aperture photometry. 

Consider multiple images ? of a single lensed quasar. Let mo{X,t) be the intrinsic 
quasar flux at time t, expressed in magnitudes at observed wavelength A. The redshifted, 
extincted flux of image i, is then 

■ A \ 



mi(A, t) = mo(A, t) - Mi{X, t) + Ei R 



zi. 



where Mi{X,t) and Ei = E{B — V) are the magnification (in magnitudes) and extinction of 
image i, and -R(A/(1 + zi)) is the extinction curve at the lens redshift zi. The magnification 
Mj(A, t) may depend on wavelength and time due to microlensing effects (Wambsganss 
2006 and references therein). By measuring the magnitude differences as a function of 
wavelength for each image pair (labeled i and j), 

mj{X,t)-mi{X,t) = AM{X,t) + AE rI J (2) 

we constrain the relative magnifications, AM(A,t) = Mj{X,t) — Mi{X,t), the extinction 
differences, AE{B — V) = Ej{B — V) — Ei{B — V), and the mean extinction curve R{X). 

For the extinction law R{X) we used either the Cardelli et al. (1989; hereafter CCM) 
parametrized models for the Galactic extinction curve or the Fitzpatrick & Massa (1990; 
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hereafter FM) model with its parameters set to match the average extinction in the SMC 
(Gordon et al. 2003). The main difference is that the Galactic models have a strong 2175A 
absorption feature while the SMC models do not. One way to confirm the presence of 
extinction is to estimate a dust redshift Zd (Jean & Surdej 1998, Falco et al. 1999), the 
redshift at which the extinction curve best fits the data, and show that it agrees with the 
observed lens redshift zi. We assume that the extinction law is the same for all images. 
Generally one image dominates the extinction and this assumption is unimportant, but it 
can be an issue if all images are significantly extincted (see Muiioz et al. 2004 and McGough 
et al. 2005). For models assuming there is only extinction, we fit the data with a single 
AM(A, t) = AM, a common differential magnification for all wavelengths which removes 
any effects from the magnifications of the macro model and most of the effects of source 
variability. 

The second physical effect in Eqn. 2 is the chromatic microlensing produced by the 
AM(A, t) term. Because of the structure in the microlensing magnification patterns, the 
changing size of the disk with wavelength changes the magnification. Since the observer, 
lens, stars and host galaxy are all in relative motion, this magnification then changes with 
time. In our present study we will examine this using simulations. Based on the properties 
of the macro models for the lens geometry, we generate magnification patterns using the 
approach of Mediavilla et al. (2006), assuming that fraction a = 0.1 of the surface density 
is in stars (i.e. that the surface density is dark matter dominated, see Kochanek et al. 2006, 
Mediavilla et al. 2009, Pooley et al. 2009, Morgan et al. 2010, Mosquera et al. 2011) and we 
simply use M = IM© stars. We then convolve the patterns with Gaussian intensity profiles 
to model the quasar accretion disk, I{R) oc exp(— i?^/2r^) where Vg oc A^ characterizes the 
disk size at wavelength A. These sizes can be rescaled to a different microlensing mass as 



M/Mq. We make many random trials fitting the data as a function of Vs and p, and then 
use Bayesian methods to estimate the size r^ and the scaling exponent p for either linear or 
logarithmic priors on r^ and linear priors on p, as explained in detail in Mediavilla et al. 
(2011). 

The last point we note is that our data are obtained at a single epoch, so our fiux 
ratios are really comparing mo(A,t) — mo(A,t + At) where At is the time delay between 
the images. This means that intrinsic source variability combined with the time delay 
between the images can lead to wavelength dependent changes in the fiux ratios which we 
will ignore by assuming that mo(A,t) — mo{X,t + At) = 0. Particularly in the estimates 
of extinction, we will see negligible effects because the parameter AM for the difference in 
the macro model magnifications also captures any achromatic effects from ignoring time 
variability. For most of the lenses we consider, these changes will be small, as can be seen 
from the empirical quasar variability models of MacLeod et al. (2010). Yonehara et al. 
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(2008) based on the ensemble SDSS quasar structure functions (Van den Berk et al. 2004; 
Ivezic et al. 2004) estimated that there would be typical shifts of ~ 0.1 mag in single epoch 
observations, but that the changes in colors would be significantly smaller because the color 
changes associated with quasar variability are far smaller than the overall variability. We 
can compensate for this problem by using modestly larger uncertainties, but it is really only 
an issue for systems with long time delays. The worst case is HE 1104—1805 which has a 
relatively long time delay of almost 6 months. If we follow the procedures of Yonehara et 
al. (2008), we estimate that the time delay can produce a bias in the shortest wavelength 
filter (F330W) of roughly 0.1 mag, with a potential color change between the F330W filter 
and the H-band of only 0.05 mag. 



3. Results 

We now consider each of the systems individually. We found extinction in 
HE 0512-3329, B 1600+434 and H 1413+117 and chromatic microlensing in HE 1104-1805. 
The remaining two systems, FBQ 0951+2635 and SBS 1520+530, did not show enough of 
a chromaticity signature to perform a deeper analysis given only a single epoch of data. 
In each of the analyses, it is necessary to determine whether the filters include any broad 
emission lines, because line and continuum flux ratios can be quite different (e.g. see 
Mediavilla el al. 2005). While both are equally altered by extinction, the broad emission 
line regions are more spatially extended and hence far less affected by microlensing (e.g. see 
Abajas et al. 2002). Here we are restricted to photometry, but by tracking the filter and 
line locations and widths we can determine the degree of contamination. Fig. 1 shows HST 
images of the 6 systems. Note that at least half of them are relatively disky, which is not 
the norm for gravitational lenses. 



3.1. HE 0512-3329 

HE 0512—3329 is a two image lensed quasar with a separation of 0'.'65 and a source 
redshift of Zs = 1.565 (Gregg et al. 2000). The lens redshift is not directly measured, but 
the presence of a damped Lymana absorber (DLA) system and associated strong metal line 
absorption systems suggests that the lens is a spiral galaxy at zi = 0.93 (Gregg et al. 2000, 
Wucknitz et al. 2003). Table 2 presents the photometry for the lens galaxy as well as the 
quasar images in the CASTLES data. While the time delays are not measured, they will be 
so short given the image separation that the single epoch flux ratios will be unaffected by 
intrinsic variability. 
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As we see in Fig. 2, the flux ratios have a steep dependence on the wavelength, and the 
slope is little changed from the earlier CASTLES results or the later results from Eliasdottir 
et al. (2006). While there are offsets between the epochs indicative of microlensing, they 
show no significant wavelength dependencies. Wucknitz et al. (2003) show that the broad 
emission line fiux ratios, which should be very little affected by microlensing, show a 
wavelength dependence consistent with these trends. After correcting for extinction using 
the broad emission line flux ratios, they also flnd an effect from microlensing. 

Figure 3 shows the result of fltting the flux ratios assuming they are due to differential 
extinction. Like Wucknitz et al. (2003), and unlike Eliasdottir et al. (2006), we identify 
a weak 2175A feature. The Galactic CCM extinction curve flts poorly, with x^ = 13 for 
4 degrees of freedom (dof), an estimated RiV) ^ 0.5 and a best flt RiV) ^ 0, a region 
where the model makes no sense. The model with the weaker feature of the mean SMC 
extinction law flts far better, with x^ = 2.8 for 5 dof and AE{B — V)= 0.06 ± 0.01. If we 
allow the parameter responsible for the stretch of the bump (c3) in the FM extinction law 
to vary, we flnd a best flt with x^ = 1.1 for 4 dof and parameter c3 = 1.7 ± 0.9 conflrming 
the marginal detection of the bump by Wucknitz et al. (2003). The interpretation of the 
flux ratios as extinction and the feature as the 2175A feature seems robust since we obtain 
a dust redshift of Zd = 0.92 ± 0.15 that is in good agreement with that of the DLA and 
metal line systems at zi = 0.93 (Gregg et al. 2000, Wucknitz et al. 2003). Although 25% 
of the CIV emission line lies in the F435W fllter, we estimate that differential microlensing 
between the line and continuum of order 0.2 mag would lead to a line-induced bias in the 
estimated continuum flux ratios of only ~ 0.01 mag. This is smaller than the photometric 
uncertainties and cannot explain the observed shift of ~ 0.1 mag. 



3.2. B 1600+434 

B 1600+434 is a two image system with a separation of 1'.'4, a source redshift of 
Zg = 1.59 and a lens redshift of zi = 0.41 (Jackson et al. 1995) where the lens is a nearly 
edge on spiral (Jaunsen & Hjorth 1997). The time delay is relatively short (~ 47 days, 
Koopmans et al. 2000, Burud et al. 2000), so single epoch flux ratios will be little affected 
by intrinsic variability. Not surprisingly, it was quickly found that the image passing 
through the disk of the galaxy suffered from extinction (Jaunsen & Hjorth 1997, Falco et 
al. 1999, Burud et al. 2000). 

Figure 4 shows the magnitude differences as a function of wavelength, where the redder 
image A is the image passing through the disk of the lens. The slope of the differences is 
little changed from the CASTLES observations, but there is an offset of approximately 
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0.2 mag. Thus, as for HE 0512—3329, the dominant effect is differential extinction with 
weaker effects due to microlensing that show no obvious wavelength dependence. Figure 5 
shows a fit to the flux ratios assuming they are due to extinction, where we have used 
the 6 cm radio flux ratio (Koopmans et al. 2000) as an extinction-free anchor for the 
ratios. The data are well fit by a CCM extinction law with R{V) = 1.5 ± 0.3 and 
AE{B — V) = 0.39 ± 0.02. The structure of the extinction law is not tightly constrained 
because the 2 175 A feature is not only bluewards of our shortest wavelength filter but also 
lies on top of the Lya line of the quasar. Dai & Kochanek (2005) estimated a gas column 
density difference between the images of ~ 3 x 10^^ cm^^ based on differences in the X-ray 
spectra of the two images. Using the extinction estimate of AE{B — V) ~ 0.1 from Falco et 
al. (1999), this implied a dust-to-gas ratio that was somewhat high. However, if we adopt 
our new estimate, we find a dust-to-gas ratio of ~ 7 x 10^^ mag~^ cm^^ that is very close 
to the typical Galactic value of 5.8 x 10^^ mag~^ cm~^ (Bohlin et al. 1978). 



3.3. H 1413+117 

H 1413+117 is a four image system with a maximum separation of ~ I'.'l and z^ = 2.55 
(Magain et al. 1988). The lens galaxy is marginally detected at H-band and its redshift is 
unknown, although Kneib et al. (1998) propose z ~ 0.9 based on the photometric redshifts 
of nearby galaxies. Figure 6 shows the magnitude differences for our HRC observations, the 
CASTLES data, the Turnshek (1997) and Chae (2001) data, and the mid-infrared (ll/xm) 
flux ratios from MacLeod et al. (2009). The small shifts between the epochs appear to 
be due to changes in the fluxes of images A and D, at levels of approximately 0.1 and 
0.05 mag. The largest wavelength dependencies correspond to images A and B and the lack 
of significant changes in the colors with time indicates that they should be attributed to 
extinction. The lack of a wavelength dependence between images D and C suggests they 
are little affected by either extinction or chromatic microlensing at these wavelengths. The 
small bump in the F435W magnitude differences including image D (see Fig. 6) is probably 
due to the contamination by the Lya emission line. In Figure 6 we also see that the flux 
ratios excluding image D {rriA — fnc, Tns — trie) extend naturally into the mid-IR as might 
be expected for extinction, while the flux rations including image D (m^ — itid-, itld — mc) 
show significant shifts (~ 0.2 mag) going from near-IR to the mid-IR. This seems more 
easily explained by microlensing, where the near-IR emission is from the accretion disk 
while the mid-IR emission is from thermal dust emission on larger scales. 

We conclude that A and B images are significantly affected by differential extinction 
from the lens galaxy. Unfortunately the lack of a candidate for the 2175 A feature combined 



with the unknown redshift of the lens galaxy makes it difficult to analyze the extinction. 
If the bump feature is present in the lens galaxy, its absence in our observations implies a 
very low lens redshift (^ 0.3), as illustrated by the example for a lens at zi = 0.25 shown 
in Fig. 7. The failure to detect the lens in the V and I-band HST observations almost 
certainly guarantees that the lens redshift cannot be so low. Thus we must conclude that 
the extinction law in this lens lacks a significant 2175A feature. 



3.4. HE 1104-1805 

HE 1104—1805 is a two image lensed quasar with a relatively large separation of ~ 3'.'2 
and a source redshift of Zg = 2.32 (Wisotzki et al. 1993) and a lens redshift of zi = 0.73 
(Lidman et al. 2000). Lehar et al. (2000) modeled the system in detail using the CASTLES 
images. The time delay is relatively long (~ 162 days, Morgan et al. 2008), but based on 
the statistics of quasar variability discussed in §2 our single epoch flux ratios should not 
be strongly biased. Falco et al. (1999) modeled the flux ratios as extinction, although the 
X-ray absorption study by Dai et al. (2006) found negligible differential absorption. In fact, 
it was also clear from the later light curves (Schechter et al. 2003, Poindexter et al. 2007) 
that there was signiflcant chromatic microlensing in this system. Indeed, as Poindexter et 
al. (2008) noted in their detailed study of microlensing the relative colors of the two images 
reversed over the period from its discovery, very different from the limited color changes 
seen in the flrst three lenses we considered. Further evidence against signiflcant extinction 
is that the mid-IR flux ratios from Poindexter et al. (2007) agree well with the emission 
line flux ratios (Wisotzki et al. 1993). Figure 8 shows the magnitude differences for images 
A and B for each ACS/HRC fllter, along with the CASTLES magnitude differences and the 
mid-IR flux ratios from Poindexter et al. (2007). We can see again the change in slope of 
the wavelength dependence between the two epochs indicating the detection of chromatic 
microlensing. 

We separately modeled the two epochs of HST observations using the procedures 
from Mediavilla et al. (2011) and Mosquera et al. (2011), as briefly outlined in §2, to 
compare the results from single epoch models to the more complex light curve modeling 
procedures used by Poindexter et al. (2008). Figure 9 shows the estimates for the scale 
radius r^ in the F336W fllter (1013A in the rest frame) and the logarithmic slope p of the 
size with wavelength, r^ oc A'', for the HRC data, the CASTLES data and the combination 
of the two assuming either a logarithmic or a linear prior on r^. In thin disk theory, where 
the disk temperature proflle is T oc R^^^^, we would expect to flnd p = 4/3. Given the 
nature of the chromatic microlensing detected in the HRC observations the uncertainties 
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are substantially greater than the ones derived from the CASTLES data, with its steeper, 
monotonic variations in the flux ratios, but the estimates (see Table 4) agree within the 
uncertainties. When we combine the two results, we find p = 1.1 ± 0.6 and r^ = 4^2 
light-days for the logarithmic prior on the size, and p = 1.0 ± 0.6 and r^ = 7 ± 4 light-days 
for the linear prior. We can compare to Poindexter et al. (2008), who used a different disk 
model and normalizing wavelength, by converting the scale lengths to the half-light radii 
Ri/2 of the models since Mortonson et al. (2005) showed that different microlensing models 
will agree on the half-light radius of the distribution. We transform our r^ at A = 3363A 
to i?i/2(A4311) = 1.18 (4311/3363)P r,(A3363) at the normalizing wavelength A = 431lA 
used by Poindexter et al. (2008), where their radius Vdisk corresponds to a half-light radius 
of Ri/2 = 2.49rA4311. In addition we rescale our microlens mass to M = O.3M0 from 
M = IMq {vg oc vM) as this is closer to the expectation for normal stellar populations 
(see Poindexter et al. 2008). Figure 10 shows that our combined results are in excellent 
agreement with those of Poindexter et al. (2008). 



3.5. FBQ 0951+2635 

FBQ 0951-1-2635 is a two image lens with an image separation of I'.'l, a source redshift 
of Zs = 1.24 (Schechter et al. 1998), and a lens redshift of zi = 0.260 ± 0.002 (Eigenbrod et 
al. 2007). The time delay is short, ~ 16 days (Jakobsson et al. 2005) and several studies 
have detected microlensing variability at the level of 0.04 mag/year (e.g., Schechter et al. 
1998, Jakobsson et al. 2005, Paraficz et al. 2006). Figure 11 shows our measurements 
of the magnitude differences between the two images along with the differences found 
by CASTLES. The magnitudes of the differences are smaller than in the previous four 
systems, and it is clear that there is little differential extinction. This agrees with the 
similar conclusion of Mosquera et al. (2011) based on ground-based narrow band imaging. 
The UV data from HST allow to us cover the wavelength range where the bump feature is 
expected given the measured redshift of the lens galaxy. Although we see a small feature in 
the F250W filter, we cannot simply attribute it to the 2175A bump because it also overlaps 
the Lya emission line of the quasar. The differences between the present data and the 
CASTLES observations indicate the presence of chromatic microlensing, but the amplitudes 
are too small for single epoch microlensing models to yield significant results. Nonetheless, 
we confirm that FBQ 0951-1-2635 is a good candidate for future measurements, reinforced 
by the fact that Morgan et al. (2010) were already able to estimate a disk size based on 
microlensing in the R-band light curves of this source. 
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3.6. SBS 1520+530 

SBS 1520+530 is a two image lens with an image separation of ~ l'/6, a source redshift 
of Zg = 1.86 and a lens redshift of zi = 0.72 (Chavushyan et al. 1997). Burud et al. 
(2002) measured a time delay for the system of ~ 130 days, and Gaynullina et al. (2005) 
and Paraficz et al. (2006) observed microlensing at a level of ~ 0.14 mag over roughly 4 
years. Morgan et al. (2010) estimated an R-band half-light radius for the disk based on 
modeling these light curves. Unfortunately our HRC observations do not show a significant 
chromaticity signal to allow us to perform a deeper analysis. Figure 12 shows the magnitude 
differences as a function of wavelength for our HRC data and the CASTLES differences 
presented in Table 2. The wavelength dependent trends in the ACS/HRC observations are 
weak, indicating that the differential extinction is very low. Interpreting the CASTLES 
data is difficult because the V-band (F555W) flux ratio is so oddly different. We have 
inspected the CASTLES data several times and have been unable to find a systematic 
problem (e.g. missed cosmic rays) that would explain the discrepancies. 



4. Discussion and Summary 

The effects of both extinction and microlensing become larger as we observe them at 
shorter wavelengths. Unfortunately, the atmosphere prevents us from observing into the 
ultraviolet from the ground, and so we generally cannot observe the rest-frame 2175A region 
to search for the characteristic feature of Galactic and LMC extinction curves or to probe 
the hot regions near the inner edges of accretion disks. Here we surveyed six gravitational 
lenses with evidence for significant wavelength dependent fiux ratios in the extinction study 
of Falco et al. (1999) from the Lband into the UV (8100A to 2200A). Two of the lenses, 
FBQ 0951+2635 and SBS 1520+530, showed changes with wavelength that were too small 
to yield interesting constraints. 

It is not surprising given the selection method that three of the lenses show significant 
evidence for differential extinction between the images. We argue for extinction dominating 
over chromatic microlensing systems based on the lack of evidence for significant time 
variability in the color, although all three systems show small changes in the fiux ratios that 
are probably due to microlensing. In the case of HE 0512—3329 we find evidence for a weak 
2175A feature from the dust in the zi = 0.93 lens. For B 1600+434 we cannot quite reach 
the wavelengths needed to quantify the presence of the 2175A feature, although a CCM 
extinction law agrees with our observations, while the lack of a lens redshift for H 1413+117 
limits our conclusions. Both systems contain significant differential extinction, and it is 
likely that the dust in B 1600+434 has the 2175A feature and that the dust in H 1413+177 
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does not. 

We clearly detect chromatic microlensing in HE 1104—1805. If we estimate the 
wavelength dependent size of the accretion disk by modeling our single epoch of data or 
the earlier CASTLES data, we find compatible results. If we combine the two single epoch 
estimates, the combined result agrees with the multi-band light curve analyses of Poindexter 
et al. (2008). Modeled as a Gaussian source exp(— r^/2r^) with r^ oc A*' and normalized at 
the observed wavelength A = 3363 A we find r^ = 4^2 (7 ± 4) and p = 1.1 ± 0.6 (1.0 ± 0.6) 
for a logarithmic (linear) prior on r^. These slopes are consistent with the expected slopes 
from standard thin disk theory (T oc R~^^^ with p=4/3), but the uncertainties are too large 
to draw a stronger conclusion. 
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Fig. 2. — Magnitude differences of HE 0512—3329 as a function of the inverse of the observed 
wavelength for the new HST observations (filled squares) along with the previous CASTLES 
observations (open triangles) and the ground-based data from Eliasdottir et al. (2006) 
(asterisks). The shaded regions correspond to the wavelength location and width of the 
most prominent quasar broad emission lines. The vertical solid line indicates the expected 
position of the 2 175 A extinction curve feature based on the estimated lens redshift. The 
horizontal error bars on the HST data indicate the widths of the filters. 
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Fig. 3. — Magnitude differences as a function of tlie inverse of tlie lens rest-frame wavelength 
for the new HST observations (filled squares). The solid line shows the best fit FM extinction 
law allowing variations in the "stretch" of the bump and the dust redshift. The dashed line 
corresponds to the best fit for a Galactic COM extinction law with i?(V^) = 0.5. 
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Fig. 4. — Magnitude differences of B 1600+434 as a function of tlie inverse of the observed 
wavelength for the new HST observations (filled squares) along with the previous CASTLES 
observations (open triangles). The format of the figure is the same as in Figure 2. 
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Fig. 5. — Magnitude differences as a function of tlie inverse of tlie lens rest-frame wavelength 
for the new HST observations (filled squares) and the 6 cm radio flux (open triangle) from 
Koopmans et al. (2000). The solid line shows the best fit for a CCM extinction law. 
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Fig. 6. — Magnitude differences of H 1413+117 as a function of the inverse of the observed 
wavelength for the new HST observations (filled squares) along with the previous CASTLES 
observations (open triangles) and the mid-infrared fiux (large open symbols: squares or 
circles when de D image is present, see text) from MacLeod et al. (2009). We have also 
shown the results from Turnshek (1997) (asterisks) and Chae (2001) (open circles). In some 
cases these points are completely hidden by our new measurements. The format of the figure 
is the same as in Figure 2. 
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Fig. 7. — ACS/HRC and mid-infrared (MacLeod et al. 2009) magnitude differences mB—fnc 
and iTiA — iTic, which are strongly indicative of extinction. As an example, the solid lines 
show fits with the same CCM extinction law for a fixed lens redshift zi = 0.25. For the more 
probable, higher lens redshift, where the location of the 2175A feature should be shifted to 
the left, successful fits require extinction laws without a strong 2175A feature. 
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Fig. 8. — Magnitude differences of HE 1104—1805 as a function of the inverse of the observed 
wavelength for the new HST observations (filled squares) along with the previous CASTLES 
observations (open triangles) and the mid-IR observations (open square) from Poindexter et 
al. (2007). The format of the figure is the same as in Figure 2. 
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Fig. 9. — Probability distributions for the size of the quasar accretion disk r^ at an observed 
wavelength of 3363A and the dependence of the size on wavelength, r^ oc A^, assuming a 
linear (left) or logarithmic (right) prior for r^, and using either the new HRC data (top), 
the older CASTLES data (middle) or the combined results (bottom). From the center, the 
contours are iso-probability density contours enclosing 15%, 47%, 68%, and 90% of the total 
probability, respectively. 
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Fig. 10. — Comparisons of the half-light radii R1/2 at A = 4311A for our combined single 
epoch models (squares) and the multi-band light curve analysis (circles) of Poindexter et al. 
(2008). Open (filled) symbols correspond to logarithmic (linear) priors on r^. The dashed 
horizontal lines represent the size estimates inferred from the black-hole mass based on the 
thin disk theory (upper) or the observed I-band flux (lower) (see Poindexter et al. 2008). In 
this composition we have shifted the mean microlens mass to 0.3 M© (see text) in order to 
better compare to Poindexter et al. (2008). 
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Fig. 11. — Magnitude differences of FBQ 0951+2635 as a function of the inverse of the 
observed wavelength for the new HST observations (filled squares) along with the previous 
CASTLES observations (open triangles). The format of the figure is the same as in Figure 2. 
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Fig. 12. — Magnitude differences of SBS 1520+530 as a function of the inverse of the observed 
wavelength for the new HST observations (filled squares) along with the previous CASTLES 
observations (open triangles). The format of the figure is the same as in Figure 2. 
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Table 1. Log of ACS/HRC Observations 



TARGET DATE-OBS 


FILTER 


EXP No. of images 


(yyyy-mm-dd) 




(sec) 




HE 0512-3329 2003-08-11 


F220W 


2136 


2 


2003-08-11 


F220W 


3x712 


2 


2003-08-11 


F250W 


408 


2 


2003-08-11 


F250W 


3x136 


2 


2003-08-11 


F330W 


129 


2 


2003-08-11 


F330W 


3x43 


2 


2003-08-11 


F435W 


36 


2 


2003-08-11 


F435W 


3x12 


2 


2003-08-11 


F555W 


27 


2 


2003-08-11 


F555W 


3x9 


2 


2003-08-11 


F625W 


24 


2 


2003-08-11 


F625W 


3x8 


2 


2003-08-11 


F814W 


6x3 


2 


2003-08-11 


F814W 


18 


2 


2003-08-11 


F814W 


3x6 


2 


FBQ 0951-1-2635 2003-10-06 


F220W 


368 


2 


2003-10-06 


F220W 


2x184 


2 


2003-10-06 


F250W 


116 


2 


2003-10-06 


F250W 


2x58 


2 


2003-10-06 


F330W 


50 


2 


2003-10-06 


F330W 


25 


2 


2003-10-07 


F330W 


25 


2 


2003-10-06 


F435W 


16 


2 


2003-10-06 


F435W 


8 


2 


2003-10-07 


F435W 


8 


2 


2003-10-06 


F555W 


12 


2 


2003-10-06 


F555W 


6 


2 


2003-10-07 


F555W 


6 


2 


2003-10-06 


F625W 


8 


2 


2003-10-06 


F625W 


4 


2 


2003-10-07 


F625W 


4 


2 


2003-10-06 


F814W 


8 


2 


2003-10-06 


F814W 


4 


2 


2003-10-07 


F814W 


4 


2 


HE 1104-1805 2003-11-05 


F250W 


2525 


2 


2003-11-05 


F250W 


842 


2 


2003-11-06 


F250W 


2x842 


2 


2003-11-05 


F330W 


303 


2 


2003-11-05 


F330W 


303 


2 


2003-11-06 


F330W 


98 


2 


2003-11-06 


F330W 


2x101 


2 


2003-11-06 


F435W 


51 


2 


2003-11-06 


F435W 


3x17 


2 


2003-11-06 


F555W 


30 


2 


2003-11-06 


F555W 


3x10 


2 


2003-11-06 


F625W 


24 


2 


2003-11-06 


F625W 


3x8 


2 


2003-11-06 


F814W 


24 


2 


2003-11-06 


F814W 


3x8 


2 


H 1413+117 2003-07-18 


F330W 


238 


4 


2003-07-18 


F330W 


2x119 


4 


2003-07-18 


F435W 


30 


4 


2003-07-18 


F435W 


2x15 


4 


2003-07-18 


F555W 


18 


4 


2003-07-18 


F555W 


2x9 


4 


2003-07-18 


F625W 


16 


4 
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Table 1 — Continued 



TARGET 


DATE-OBS 


FILTER 


EXP No. of images 




(yyyy-mm-dd) 




(sec) 




SBS 1520+530 


2004-06-15 


F250W 


2484 


2 




2004-06-15 


F250W 


3x828 


2 




2004-06-15 


F330W 


561 


2 




2004-06-15 


F330W 


3x187 


2 




2004-06-15 


F435W 


141 


2 




2004-06-15 


F435W 


3x47 


2 




2004-06-15 


F555W 


90 


2 




2004-06-15 


F555W 


3x30 


2 




2004-06-15 


F625W 


68 


2 




2004-06-15 


F625W 


3x26 


2 




2004-06-15 


F814W 


62 


2 




2004-06-15 


F814W 


3x24 


2 


B 1600+434 


2003-08-17 


F330W 


4x1080 


2 




2003-08-17 


F435W 


4x498 


2 




2003-08-17 


F555W 


4x402 


2 




2003-08-17 


F625W 


4x312 


2 




2003-08-17 


F814W 


4x293 


2 
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Table 2. CASTLES Photometry 



Lens 


Component 


AR.A. (") 


ADec {") 


F555Wt 


F814W 


F160W 


HE 0512-3329 




mage A 


0.182 + 0.003 


0.621 + 0.003 


18.15 + 0.06 


16.81 + 0.08 


15.81 + 0.02 






mage B 








18.40 + 0.09 


17.28 + 0.07 


16.38 + 0.03 






lens G 


0.09 + 0.07 


0.37 + 0.10 


22.1 + 0.6 


20.9 + 0.7 


19.1 + 0.8 


FBQ 0951+2635 




mage A 








17.29 + 0.06 


16.70 + 0.03 


15.62 + 0.03 






mage B 


0.900 + 0.003 


-0.635 + 0.003 


18.32 + 0.06 


17.89 + 0.02 


16.99 + 0.03 






lens G 


0.760 + 0.003 


-0.455 + 0.003 


21.02 + 0.04 


19.67 + 0.03 


17.86 + 0.14 


HE 1104-1805 




mage A 








16.92 + 0.06 


16.40 + 0.03 


15.91 + 0.01 






mage B 


2.901 + 0.003 


-1.332 + 0.003 


18.70 + 0.08 


17.95 + 0.04 


17.35 + 0.03 






lens G 


0.965 + 0.003 


-0.500 + 0.003 


23.26 + 0.30 


20.01 + 0.10 


17.52 + 0.09 


H 1413+117 




mage A 








18.00 + 0.01 


17.77 + 0.01 


15.83 + 0.04 






mage B 


0.744 + 0.003 


0.168 + 0.003 


18.07 + 0.01 


17.84 + 0.01 


15.92 + 0.03 






mage 


-0.492 + 0.003 


0.713 + 0.003 


18.27 + 0.01 


18.06 + 0.01 


16.18 + 0.02 






mage D 


0.354 + 0.003 


1.040 + 0.003 


18.32 + 0.01 


18.15 + 0.01 


16.43 + 0.03 






lens G 


0.142 + 0.003 


0.561+0.003 


- 


- 


18.61 + 0.03 


SBS 1520+530 




mage A 








18.83 + 0.05 


17.97 + 0.03 


17.58 + 0.02 






mage B 


1.429 + 0.003 


-0.652 + 0.003 


19.29 + 0.24 


18.99 + 0.07 


18.41 + 0.03 






lens G 


1.141 + 0.003 


-0.395 + 0.003 


23.40 + 2.00 


20.16 + 0.11 


18.22 + 0.05 


B 1600+434 




mage A 








23.61 + 0.12 


21.92 + 0.10 


20.66 + 0.03 




image B 


-0.720 + 0.003 


1.183 + 0.004 


22.32 + 0.09 


21.39 + 0.03 


20.47 + 0.03 




lens G 


-0.110 + 0.003 


0.369 + 0.004 


- 


20.78 + 0.06 


18.30 + 0.13 





Note. — t For the system H 1413+117 it corresponds to the filter F702W 



Table 3. ACS/HRC Photometry 



Lens 


Image 


F220W 


F250W 


F330W 


F435W 


F555W 


F625W 


F814W 


HE 0512-3329 


A 


18.96+0.11 


18.07+0.23 


17.67+0.13 


18.67+0.03 


18.10+0.05 


17.60+0.05 


16.98+0.03 




B 


18.33+0.04 


17.74+0.02 


17.55+0.03 


18.66+0.02 


18.25+0.04 


17.88+0.03 


17.36+0.03 


FBQ 0951+2635 


A 


16.72+0.01 


16.36+0.03 


16.60+0.02 


17.80+0.04 


17.48+0.03 


17.14+0.03 


16.82+0.03 




B 


17.97+0.03 


17.70+0.08 


17.82+0.01 


19.00+0.08 


18.71+0.04 


18.40+0.10 


18.12+0.05 


HE 1104-1805 


A 








17.25+0.05 


17.81+0.07 


17.57+0.10 


17.33+0.06 


16.85+0.05 




B 








18.20+0.09 


19.04+0.12 


18.90+0.11 


18.71+0.12 


18.18+0.07 


H 1413+117 


A 






20.84+0.03 


18.97+0.07 


18.61+0.06 


18.20+0.09 


17.75+0.02 


17.70+0.03 




B 






21.40+0.12 


19.36+0.08 


18.91+0.12 


18.48+0.07 


17.95+0.05 


17.88+0.01 




C 






20.45+0.01 


19.00+0.01 


18.84+0.03 


18.53+0.04 


18.13+0.05 


18.10+0.02 




D 






20.78+0.10 


19.40+0.06 


19.20+0.09 


18.69+0.02 


18.26+0.04 


18.25+0.01 


SBS 1520+530 


A 






18.23+0.07 


17.86+0.03 


18.94+0.02 


18.73+0.03 


18.52+0.01 


18.12+0.04 




B 






19.40+0.10 


18.93+0.03 


19.94+0.02 


19.66+0.04 


19.46+0.01 


19.10+0.04 


B 1600+434 


A 








25.68+0.47 


25.36+0.11 


24.63+0.10 


23.67+0.02 


22.68+0.03 




B 








22.55+0.17 


23.49+0.25 


23.07+0.11 


22.44+0.03 


21.76+0.06 
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Table 4. Quasar Accretion Disk Measurements for HE 1104—1805 





ACS 


CASTLES 


ACS 


X CASTLES 


Logarithmic prior 


rs (light-days) 
P 


6tl 
1.8±0.8 


l.OiO.7 




1.1±0.6 


Linear prior 


rs (hght-days) 
P 


12±6 
1.8±0.8 


7±4 
0.9±0.6 




7±4 
1.0±0.6 





Note. — Tg is the size of the quasar accretion disk modeled as 
a Gaussian {I{R) oc exp(— i?^/2r^)) at the observed wavelength 
A = 3363 A and p is the power law of the size variation with 
wavelength (^^(A) oc A^). 



